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Abstract 

The explicit form of the interface equation of motion derived assuming a minimal surface 
is extended to general bicontinuous interfaces that appear in the diffusion limited stage of the 
phase separation process of binary mixtures. The derivation is based on a formal solution of 
the equivalent simple layer for the Dirichlet problem of the Laplace equation with an arbitrary 
boundary surface. It is shown that the assumption of a minimal surface used in the previous 
linear theory is not necessary, but its bicontinuous nature is the essential condition required 
for us to rederive the explicit form of the simple layer. The derived curvature flow equation 
has a phenomenological cut-off length, i.e., an 'electro-static' screening length. That is related 
to the well-known scaling length characterizing the spatial pattern size of a homogeneously 
growing bicontinuous phase. The corresponding equation of the level function in this scheme 
is given in a one-parameter form also. 

1 Introduction 

The ordering process associated with the first order phase transition has been investigated for a 
long time since the early works by van der Waals appeared at the end of the last century. [|T| It has 
occupied a major position together with the dynamical critical phenomena as a subject not only 
of the non-equilibrium statistical physics but also of the field theory. A challenging feature of it is 
the strong nonlinearity which causes a variety of evolving spatial patterns. This has initiated the 
interesting concept of 'pattern formation' in far-from-equilibrium thermodynamic systems. 

Complicated patterns, at a glance, are seen even in the most simplified theoretical models of 
a binary mixture, e.g., in the time-dependent Ginzburg-Landau (TDGL) model for the conserved 
order parameter (COP), or the Cahn-Hilliard (CH) model. The nonlinearity in these systems 
finally gives rise to a spatial singularity of the order parameter gap, i.e., the interface with a sharp 
profile that separates distinctly the coexisting couple of phases. This is a primary example of the 
so-called 'topological defect' appearing in continuous fields. The main subject for us in this stage is 
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to reduce a degenerate evolution equation for the interface from the basic TDGL or CH equations 
for the bulk field. 

In spite of remarkable success of the intuitive, droplet theories on a dilute mixture, i.e., some 
pioneering works in the middle of this century, f8fl-|[10|| the interface evolution equation for more 
general systems has not been obtained yet in an explicit form. The most general form we have 



found so far is an integral equation obtained by Kawasaki and Ohta. [11, 12| This is a kind of the 



curvature flow equation which represents the interface velocity with the mean curvature implicitly in 
contrast to the Allen-Cahn equation [|l!|] for a system of non-conserved order parameter (NCOP). It 
has been shown |]14| that this is the equivalent simple layer equation |15] for the Dirichlet problem of 
the Laplace equation that is derived as a quasi-static approximation for a Stefan problem, |f[6f 



i.e., 



a diffusion problem of an autonomous boundary condition evolving through the diffusion process 
itself. For some simplified systems, e.g., a planar interface or a thin system of spherical droplets, 
it is straightforward to derive the curvature flow equations, if they are considered as electro-static 
potential problems. The explicit potential solutions for such systems can be found in every ele- 
mentary textbook on electro-statics. Jl7j On the contrary, we will soon come up against a serious 
difficulty of a complicated geometry of the boundary condition when we try to find out explicit 
potential solutions for more general systems. For example, in a nearly symmetric mixture there 
appears a random bicontinuous structure, i.e., the so-called sponge phase. This has an infinitely 
multiply-connected topology and is characteristic of a three dimensional system. That contrasts to 
the case of a lamellar phase, which is essentially equivalent to a planar system. 

In the previous work[14[] the author assumed an almost minimal surface for such a system mainly 



by intuition gained from some simulated pictures [|T^] - ]20|| , which just remind us a periodic minimal 



surface such as the Schwarz lattice structure, PT] at least in a local view. The theoretical foundation 
for this assumption was that the mean curvature is included linearly in the boundary value itself 
as the Gibbs-Thomson condition already, and then it can be neglected within the limits of a linear 
theory with respect to the small mean curvature. This situation is similar to that of the decay 
process of a perturbative deformation on a planar interface. We used the common feature of zero 
mean curvature of the unperturbed basis in both systems. An important difference is that the 
sponge phase has an evident characteristic length, which was used effectively as the mean electro- 
static screening length in the previous work. 

In the present paper it is shown that the assumption of a minimal surface on the unperturbed 
basis is not necessary in deriving an explicit solution from the integral equation and the same one- 
parameter curvature flow equation is rederived for more general bicontinuous systems. We need 
the bicontinuous nature of the interface only. Here 'bicontinuous' means that the interface divides 
the whole R 3 space into a couple of ( not necessarily symmetric) subspaces connected in each like a 
planar interface does, and not into subspaces more than three. This is the only condition required 
to make up a formal expression for the equivalent simple layer in terms of the potential problem. 

In §2 the integral equation for the interface velocity is rederived on the basis of Onsager's 
variational principle. It is shown that the principle of minimum dissipation is an example of the 
so-called gradient dynamics and is very useful for such an actual purpose to reduce the degenerate 
interface equation of motion correctly from the basic transport equation for the bulk field. In §3 a 
formal explicit expression for the simple layer on an arbitrary connected surface is derived referring 
to a solvable problem of a plane boundary. The equation of the level function that corresponds to 
the present curvature flow equation is discussed in §4. 
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2 The Onsager Principle and Basic Equations 



In the original work by Kawasaki and Ohta[TT|, [12] a kind of the path- integral method and a 
variational principle associated with it were used. It seems instructive for us to follow it in terms 
of the familiar Onsager principle of minimum dissipation [[22] with reference to the recent gradient 
dynamics. To begin with, let us survey the Onsager variational principle. 

Let x = {x-i} be a set of generalized thermodynamic variables, their phenomeno logical transport 
equations being given by 

Xi = ^2 LijXj, (1) 

3 

where {Lij} are the Onsager coefficients and X = {Xi} is a set of generalized thermodynamic 
forces, which are defined by using the entropy function S as 

dS 

X = \7 X S or X t = — . (2) 

Onsager proposed a minimum principle on the basis of the symmetric property of {L^}, i.e., = 
Lji, and its positive definiteness as follows: Define two kinds of dissipation function and a Lagrangian 
by 

&[x,x] = 7j Yl L^ijX&j, (3) 

Z i,3 

9[X,X] = ^L ij X i X J , (4) 



and 



C[x,X] = $[x,x]-S[x,X] + V[X,X] 

= i^L^ijixi -J2 L ikXk)(ij -J2 L H X i)i ( 5 ) 

i,j k I 



where 



S[x,X] =Y,XiXi. (6) 

i 

Then the phenomenological equation Eq.(|l]) is given by 

5C = (minimum) with respect to x. (7) 

This variational principle was confirmed later as a problem of the most probable path in the path 
integral method of the linear fluctuation theory |23|, [24]] It has been understood widely that the 
principle itself is nothing but a formal theory as is suggested in the above survey and is no use for 
an actual purpose to find out the transport equation itself for a given system. 

Recently, mathematicians have introduced a notion of 'gradient dynamics'. [25] In this sense, the 
above Onsager's phenomenology is a Lagrange multiplier version[] of a gradient dynamics defined 
by an inner product, 

O, v) =J2 L^ijUiVj = 2$ [it, v}. (8) 

1 Th.e original plan is that a dynamics is to be defined by maximizing S = x - VxS with a constraint of a properly 
chosen inner product (A, A) = constant. 
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Mathematicians may define an arbitrary dynamics by using an arbitrary inner product. That gives 
us a good hint: That is, it can be reasonably expected that we may find out the correct interface 
dynamics easily, when we rewrite the inner product, i.e., the physical dissipation function for the 
bulk into an interfacial form. 

Let {s(r)} be the scalar field whose evolution is described by the TDGL equation of COP type, 
or equivalently the CH equation, 

!•<-•'> =- iv, j£j *«*■»>■ < 9 > 

In this case S'({s(r)}) is not the entropy but related to the free energy functional F({s(r)}) by 
S = —F/kT, where T is the temperature and k the Boltzmann constant and 

F({s(r)}) = J dr^-{- l -s{rf + -^A r f + (10) 

The equation (||) has a uniform equilibrium solution, 

s(r) = ±so, (11) 

and a planar interface one, i.e., the so-called kink solution, 

z 

sk(z) = s tanh— , (12) 

where z is a normal coordinate perpendicular to the interface. Here £ corresponds to the thickness 
of the interface. The parameter xo is related to the susceptibility by 



_! _ S 2 F 
^° 5s(r) 2 



±s 



(13) 



From Eq.(|P2"D we get an expression of the surface tension, |]J 



e r°° 2 e 2g 

a = — / (Vs K ) dz = — , (14) 

as the excess free energy stored in the interface layer where s ^ ±Sq. 

Now let us derive the interfacial version from the bulk dissipation function defined by 

where Gq(v — r') is the inverse of the Onsager coefficient —V 2 in Eq.(||) defined by 

V 2 G (r-r') = -5(r-r'), (16) 

which is the Green function, i.e., the Newton (or the Coulomb) potential. If we assume that the every 
element dA of the interface S at a propagates without deforming its profile, it is straightforward to 
show that 

^P-dz = -2s v n (a). (17) 

3 Ot 
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from a geometrical consideration, where v n (a) is the interface normal velocity. Thus, the interface 
version of the dissipation function Eq . fljjD is given by 



* (&o) s 



)? / / G {a-a')v n {a)v n {a')dAdA', (18) 
2L JsJs 

where the volume elements dr and dr' in Eq.flTo]) are replaced by dzdA and dz'dA'. On the other 
hand, by definition of the surface free energy a, the interface version of S is written as 

S = ~ <f aH(a)v n (a)dA, (19) 

according to the formula on the change of the surface area caused by the normal displacement 5z 
of the surface S, 

5 j> dA = j> H{a)8zdA, (20) 

where H(a) is the mean curvature of the surface at a. Then, by using the variational principle 
<5($ — S) = with respect to v n (a) we obtain an integral equation, 

/ G (a - a')v n (a')dA' = [H(a) - H], (21) 

Js 6 

which is just the equation obtained by Kawasaki and Ohta. Here a new notation, 

D = Tfr (22) 



is introduced for the simplicity together with the relation Eq. (|T4]) . Note that a constant, H, is 
incorporated as the Lagrange unknown multiplier in the present variational principle, corresponding 
to the constraint of order parameter conservation, 

^v n {a)dA = 0. (23) 

If the inverse r (a, a') of the integral kernel Go(a — a') satisfying 

j> r (o, a")G (a" - a')dA" = S(a - a'), (24) 

were given, the constant H would be determined by a weighted average, 

H=ii r (o, a')H(a)dAdA'/ 1 1 r (o, a')dMA'. (25) 

Apparently the integral equation Eq . (|2T|) having a kernel of the Newton potential is the same as 
that of the equivalent simple layer for the Dirichlet problem of the Laplace equation. Let us relate 
it to the diffusion equation in the followings. 

Because the order parameter deviation, Ss, from the saturated level ±s at each phase is expected 
to be very small in the region out of the boundary layer (\z\ > £), the chemical potential can be 
approximated by 

d 2 F 

li ~ —y Ss = Xo Ss. (26) 

OS ±s o 
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Therefore, the limiting process in this region must be an ordinary diffusion process defined by 

= DV 2 ^j, (D = L/kTxo) (27) 

where a normalized variable, ip = Ss/sq, is introduced. By using a kind of the boundary layer 
method we find a boundary condition, 

2s /i(0) = cxH(a), (28) 

that is, the well-known Gibbs-Thomson condition for the curved interface having the mean curvature 
H = — V ■ n, where n is the normal unit vector of the interface. Then, the boundary value for ip 
must be given by 

^s(a) = \iH{a). (29) 
It should be noted that the boundary S deforms itself via this diffusion with a normal velocity, 

v n (a) = ^A[n-V^] a , (30) 

where A[...] a denotes a gap of the enclosed quantity through the boundary layer at a. These 
equations make a Stefan problem with a time- dependent, autonomous boundary condition. 

In the late stage of the phase separation process where the characteristic length, A, of the spatial 
pattern is sufficiently greater than the interface thickness £, the propagating velocity f n (~ DVip ~ 
D£H/\ ~ D£/\ 2 ) and the diffusion velocity Ud(~ DVip/ip ~ D/X) satisfy the condition 

^ ~ | « 1. (31) 
v D X 

Therefore, a quasi-static approximation, 

V 2 ^ = 0, (32) 

is applicable in this stage. Thus, the above Stefan problem becomes a Dirichlet problem of the 
Laplace equation. The propagating velocity Eq. ([30) together with the boundary condition Eq. ( |2"9"| ) 
is now given by the gap of the 'electric field' across the boundary, i.e., the equivalent 'simple charge 
layer', which should satisfy the integral equation Eq.([H|). The additional parameter H should be 
regarded as a compatibility condition in this framework as follows: When {ip{r)} is the solution 
for the boundary condition {-0g(a)}, {ip{r) + c} is the solution for another boundary condition 
{ifjs(a)+c} with an arbitrary additional constant c. Obviously the family of the boundary conditions 
{ip s (a) + c} have the same value of the field gap A[n • Vip]a or the simple layer solution. If the 
equivalent simple layer solution for this family exists, it determines a unique boundary value given 
by the left hand side of the integral equation. Then the right hand side should have an adjustable 
constant for the compatibility when an arbitrary boundary condition {tps(a)} is given. Note that 
this is true only if no flux lines escape out of the system into infinite points as is mentioned below. 
That corresponds to the conservation condition Eq.(^), because of the Gauss theorem, i.e., the 
total escaping flux being given by the total charge, 



j> v n (a)dA, 
J s 



(33) 
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included in the system. Thus the constant H must be determined by the conservation condition 
again. An exceptional example not covered by this rule is a problem of an isolated system composed 
of finite closed surface (s), where the flux may escape out. Of course we know that this is an ordinary 
case in the usual electro-static problem. For the simplicity suppose a finite sphere of radius R and 
a constant boundary value on it- The simple layer in this case is given by A[n • Vip] = ips/R, 
because ij)(r) = Rips/ r f° r r > R and ip(r) = ips for r < R. Evidently, another boundary condition, 
ips + c ) has a different value of the simple layer, (ip s + c)/R. Thus, for the problem of isolated closed 
surface(s), the compatibility constant is not necessary. In contrast to it, for the extraordinary case 
of infinitely extended systems such as a system of spheres scattered homogeneously in the whole R 3 
space or a bicontinuous system spreading over the whole R 3 space, etc., we need the compatibility 
condition. In these systems the constant boundary values induce a unique constant potential in 
both sides of the interface(s). 



3 Derivation of the Explicit Solution 

3.1 A planar boundary 

First let us consider a Dirichlet problem of the Laplace equation, 

V 2 ^ = 0, (34) 

when the values {ips(x, y)} on the xy-plane, S, are given. Let dufi',!/'^) be the surface charge 
density on the grounded conductor S induced by a probe charge —1 located at a referred point 
(0, 0, z). Then the potential ip at the referred point (0, 0, z) is expressed in an explicit form as 

ip(0,0,z) = / / a ind (x',y';z)i! S (x',y')dx'dy', (35) 



by using the elementary method of the Green function. The induced charge density <r in d on the 
plane boundary can be found in every textbook on electro-statics as the simplest example of the 
mirror image method. By using the actual expression for it together with the Taylor expansion, 

rfatf, y') = exp(a' ■ V s )^s(0, 0), (36) 
where a' = (x',y f ) and Vs = (d/dx,d/dy), Eq.([35l) is rewritten as 



^(OAz) = J J 2<x , 2 - yl2 - z2)3/2 exp(a' ■ V s )^s(0, 0)dx'dy' 

roo 2,0! 

= / da ' t al 2W2 J o(o^Vsl) Vs(0,0) (37) 
Jo [a z + z z ) 6 ' z 

= exp(-2;|zV s |) Vs(0,0), 

where J n (r) is the usual Bessel function of n-th order. Then, by differentiating it with respect to 
z, we find 

n-V^U +0 = -KV s |^s(a)- (38) 
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Note that this result itself means merely a transform of the boundary condition from a Dirichlet 
type into a Neumann type. Combining with the same, symmetric result for the opposite side, 
z — ► —0, we obtain the formula for the equivalent simple layer, 



(T, 



eq 



» = -A[n-V^] a = 2|zV s | [Vs(a)-$s], (39) 



by means of which we can write the potential in another explicit form, 

^(r) = I G (r-a')a cq (a r )dA'. (40) 

w S 

The compatibility constant mentioned in the last part in §2 is incorporated here, although it would 
disappear after the differentiation in the present case. It should be noted that, from the view-point 
of the Green function method this term is required, in principle, for the condition ip — > when 
\z\ — ► oo on applying the Green theorem. 



The fictitious operator |zVs| was introduced by Ohta and Nozaki[|26| in a perturbational theory 
for a planar interface. They defined it by a Gaussian integral. However, in the context of the 
present derivation it should be interpreted as the following limit, 

|zV s | = lim V(-A 2 V S 2 ), (41) 

A— >oo ,\ 

where Vs 2 (= d 2 /dx 2 + d 2 /dy 2 ) is the surface Laplacian and the function V(Q 2 ) is defined by 

CO 1 

= HQ) - QUQ) + Q [ Q MQ')dQ'- ( 42 ) 

Jo 

The asymptotic behavior for large \Q\ is evaluated as 

V{Q 2 ) ~ |Q|, (43) 

which results in the well-known |g| 3 dispersion relation|27|]-|29| for the relaxation mode on a planar 
interface obeying the decay law ~ exp(— D£q 3 t/3). The formula Eq. (f42"D was derived in the previous 
paper by introducing a cut-off length or the upper limit A for the integration in Eq. (p?|) . This 
time the cut-off is temporarily required to integrate the series expansion for the Bessel function 
Jo(a'|iVs|), or for exp(a' ■ Vs), term by term. It will be regarded as a physical screening length and 
has an essential role in the final formula for bicontinuous systems in the followings. 

The above formulation is applicable to the two dimensional system also, where the function 
V{Q 2 ) should be replaced by 

v(Q 2 ) = --t T o — W^"^ n 

7r^ (2n-l)(2n)! 

2 rQ sin O' 

= ±[cosQ + Q S -^-dQ'] , (d = 2). (44) 

7T JO Q 

The asymptotic behavior is the same as Eq . (f43l) . Here it should be noted that the operator \id/dx\ 
in this case is neither d/dx nor \df /dx\ if operated on f(x). 
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3.2 General curved surface boundaries 



Now let us discuss the general case of a parametric, curved surface S. Let (ui, u 2 , u 3 ) be a curvilinear 
orthogonal coordinate system, where the surface S is defined by u 3 = 0. The Laplacian in this system 
is given by 

1 3 d 'y 8 
7 dui g { 2 dut 

where 

g i = \ 7 ^\ = l/\Vu l \, (46) 

OUi 

is the linear metric (i-th element of the diagonalized metric tensor) and 7 is the Jacobian defined 
by 7 = 9x929?,- We need the infinitesimal vicinity of S only in the followings. Then it is always 
possible without loosing the generality to assume g 3 {v,\,u 2 ,u 3 ) = 1 everywhere by introducing a 
family of parallels |30| of S for the simplicity. In addition, let (ui,u 2 ) be the orthogonal coordinate 



associated with the directions of principal curvature on S. Suppose that a probe charge —1 is located 
at a referred point (0, 0, z) in this curvilinear coordinate and the surface S is grounded. The Green 
function G(ui, u 2 , u 3 ; z) of the Laplace equation for this boundary condition obeys the Poisson 
equation, 

£ G = Kui)5(u 2 )5(u 3 - z), (47) 



:1 du i dui 



with the constraint, 



G(ui,u 2 ,0;z) = 0. (48) 

Here new parameters e\ = g 2 /gi,e 2 = 91/92 arid e 3 = g%g 2 are introduced for convenience. Note 
that the denominator 7 in Eq . (f45|) canceled out with the Jacobian, which had appeared in the right 
hand side of Eq. (T4"T|) . This can be regarded as an electro-static problem with a planar boundary 
condition in a Euclidean space (ui,u 2 ,u 3 ) having an anisotropic, heterogeneous dielectric tensor, 
if the whole parametric space (ui,u 2 ,u 3 ) corresponds to the original physical space, and vice versa. 
For the time being let us assume it on condition that we need only the limit z — > 0, and let it be 
discussed later. 

Therefore, the elementary mirror image method is applicable to this case, taking account of the 
refraction of the flux lines, 

_ / dG dG 9G\ ,. n . 
D= (-e 1 — ,-e 2 — ,-e 3 — ). (49) 



du\ du 2 du 3 > 

For the present purpose to find n-Vip on S we need only the flux lines that run along the infinitesimal 
vicinity of S in order to calculate the induced charge <7i n d, because the fields produced by the probe 
charge at z and by its image at — z can be calculated separately. 

The method of a reference frame used in the previous paper to calculate the refracted flux 
assuming a minimal surface can be extended to more general surfaces as follows: If we take a 
normalization, ej(0, 0, z) = 1, for the simplicity, the flux that arrives in the area element du\du 2 on 
the plane S is related to the solid angle dVt when the flux started from the source point (0, 0,z), 
that is, D 3 du\du 2 = dfl/ATr. Let du\du 2 be the area element that dfl would cut from the plane S 
when it were extended straight without refraction. Then the solid angle dfl is given by 

dQ = t~2^~2 , 2 ^ l2 duidu 2 . (50) 
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Thus, the required surface charge induced in the area element dA on the original curved surface S 
is finally given by 

dG z 
o"ind(wi, u 2 ] z)dA = -- — dA = 2 - 2 — ^Tfyduidu 2 , (51) 

<9u 3 «3=0 27r(Mi + u 2 + z 2 ) 3 ' 2 

where the relation dA = e^du\du 2 is used and a factor 2 caused by the image charge is introduced. 
Note that the actual expression for the refraction law of the flux, or that for the mapping (ux, u 2 ) — > 
(u 1 ,u 2 ) are not necessary in the above calculations. 

This is the method of the pseudo-conformal transformation introduced in the previous paper. 
Thus the assumption of a minimal surface used there has none of special meaning for us no longer. 

3.3 The bicontinuous phase 

The remaining procedure is almost the same as that for the planar boundary except for the following 
three facts: First, it should be noted that the operator we will face in the expected formula related 
to Eq.(^) after the Taylor expansion must be 

d 2 + d 2 
du\ du 2 

However, in the limit z — > 0, this can be replaced by the desired surface Laplacian defined by 

g t g 2 v dux 9i dux du 2 g 2 du 2 ' 

as is shown in Appendix A. Note that this fact does not mean a conformal mapping as is mentioned 
there. 

Second, there still remains the question whether the whole physical space is covered by a para- 
metric space (ui,u 2 ,us) or not. Apparently the answer is no in general except for simple curved 
surfaces having the same topology as that of a plane. In a complicated sponge structure the region 
described by a given set of parameters (ux,u 2 , U3) must be limited in the local cave around the 
referred point. In order to avoid this difficulty, let us introduce an upper-limit A for the integration 
in Eq. fl37|) . This cut-off length was defined intuitively as an electro- static screening length, i.e., 
the effective diameter of local caves of the sponge-shaped conductor in the previous paper. In fact 
we have such a unique (but time-dependent), well-defined characteristic length in the bicontinuous 
phase that grows up starting from a quenched, homogeneous mixtures. It should be noted that we 
have another difficulty that at least one of the parameters (ui,u 2 ) may become multi-valued, for 
example, on a part of the surface having a rotational symmetry like a catenoid. It seems that we 
have no problem in applying the image method to this case if we do not introduce any cut-off, as 
is exemplified using a two-dimensional solvable problem of a circle boundary in Appendix B. The 
cut-off procedure, however, causes literally cut-off of the multi-valued part in this case. Fortunately, 
in addition to the fact that we need only the very vicinity of S, the weight of the contribution of 
the boundary values at the points away from the referred point decreases as <7i n d ~ 1/f 3 , where 
f = [ui + u 2 ) x l 2 is essentially the distance along S as is suggested by the equality between the 
operator Eq. (|52|) and Vs 2 . Thus, the cut-off procedure may be introduced without ruining the 
approximation. 
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The last point is the bicontinuous nature of the interface S. This is related to the equivalence of 
the boundary condition between the inner and the outer problems of S. For example, let us consider 
a non-bicontinuous system composed of several closed surfaces, Si,S2, ... • Evidently, each inner 
problem associated with each Sj has a single boundary Sj itself. On the other hand the sole outer 
problem has a united boundary Si +S2 + ••■ • Thus the boundary conditions are not equivalent to 
each other. In contrast to it, we have just two regions separated by a single sheet of the interface S in 
the bicontinuous phase that is extended infinitely or is connected periodically in all three directions, 
and the difference between the 'inner' and 'outer' problems has no sense there. Both problems have 
exactly the same boundary in this case. Then the induced charge densities for both problems must 
coincide with each other when z — > ±0 as is seen in the above derivation, even if the both bulk 
regions are asymmetric at an average. Therefore, the values of the potential gradient n • Vip(a) 
constructed with the common boundary values {ips(a)} on S become equivalent one another. 

Thanks to these facts we can construct the equivalent simple layer by using Eq . fl39|) . Thus the 
explicit expression for the curvature flow equation, 



is obtained again for general bicontinuous phases. The parameter A(t) is the time-dependent char- 
acteristic length, which is to be determined self-consistently by this equation. On the assumption of 
scaling, i.e., the time-dependent similarity law, the well-known time-dependence A(t) ~ t 1 ^ 3 may be 
found by using a kind of dimension analysis on this equation. Note that the compatibility constant 
H in the right hand side is given by the simple surface average in this approximation, i.e., 



The first term expresses the mean-field evaporation-condensation process (~ {H{a) — H)), due 
to the imbalance of the mean curvature around its average value H. The second term may be 
interpreted as a kind of surface diffusion (~ —X7sH(a)) due to the local nonuniformity of the 
curvature. On the contrary the asymptotic value of V{Q 2 ) for Q — > ±00 is given by Eq . fl43|) , i.e., 
V(Q 2 ) ~ \Q\, which corresponds to the |g| 3 relaxation mode (~ — |Vs|i?(a)). Then the short wave- 
length local fluctuations (|g|A(t) 3> 1) on the interface must decay rapidly and we can expect that 
the interface is always smooth within a spatial scope of X(t) almost everywhere. 

The above behaviors of V(Q 2 ) are shown in Fig.l. The asymptotic estimation (~ \Q\) is practi- 
cally applicable for Q > 2. However, the important modes for the interface evolution are included in 
the small wave-number region, |gA(t)| < 1, where V(Q 2 ) deviates significantly from its asymptote. 

4 Discussions 

The Onsager principle of minimum dissipation seems to be really useful in such a practical problem 
to reduce the interface evolution equation from the bulk equation, at least if the effect of thermal 




(54) 




(55) 



For small Q, the function V(Q 2 ) is expanded asBT] 




(56) 
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fluctuations can be neglected. In the unstable phase separation process the main role of thermal 
fluctuations is creating nuclei or microscopic droplets at the very early stage. This effect may be 
incorporated as the initial condition in the picture of the interface dynamics. Another example is 
the NCOP system. The inverse of the Onsager coefficients in this case is merely L~ 1 5(r — r'). Then 
the dissipation function is given by 

*--rJ^--mi^ afdA < (57) 

where the definition of the surface tension Eq . fll4|) is used. On the other hand, S becomes 

S = ~ <f [oH(a) - 2s h}v n (a)dA, (58) 

where h is the uniform external field that should be included in the free energy in the following 
form, 

— J hs(r)dr, 

for a NCOP system. Then using the variational method we find 

v n (a) = -^-[H{a) - ^ ]. (59) 

That is the well-known Allen-Cahn equation. In addition, putting H = (d—l)/R for the droplet 
system, we obtain a radius evolution equation, 

. (d - i)Le 1 1 i \ , Kn , 

where R c — (d — l)a/2soh is the radius of the critical droplet. Thus we can get these important 
formulae straightforwardly and correctly by using the Onsager principle. 
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The one-parameter curvature flow equation for general COP systems is rederived without as- 
suming a minimal surface. The bicontinuous nature is shown to be indispensable in the present 
derivation. Of course the result is applicable to a planar interface if we adopt an infinite cut-off 
length, A —>■ oo. Especially it should be noted that it gives an exact explicit solution of the Dirichlet 
problem for an arbitrary curved surface having the same topology as that of a plane. Such formulae 
have not been found in the textbook on the electro-statics or on the potential problem. 

On the contrary, this result cannot be applied to a system of scattered spheres because it has 
not the bicontinuous nature, although it has been one of the simplest examples to study. Of course 



a more accurate argument on this system f32fl is possible by applying the standard electro-statics 



directly to Eq.fl21~D as was performed by Kawasaki and Ohta. |L2| Suppose each sphere of radius R 



and charged by Qi is sufficiently separated from each other. Then Eq . (|2T|) for the potential on the 
surface of the sphere i becomes 



AirRi f^^irRij 6 R, 

where Rij is the distance between a couple of spheres i and j and is assumed as Rij 3> Ri, Rj. Then 
the surface charge density given by Ri = Qi/AnRi 2 satisfies 

R +J2 Rj2 R - D ^ 1 [ 1 1 ] (62) 
. ,j RiRij 3 Rt R c (t) Ri 

The second term in the left-hand side shows plainly the long-range interaction of the interface 
through the diffusion process. When this term is neglected for a sufficiently thin system, this 
becomes the Lifshitz-Slyozov equation, where the critical radius R c (t) in the right-hand side is 
to be determined by 

R c (t) = (R)= y £Ri/T, 1 > ( 63 ) 

according to the conservation condition YU ^Ri'Ri = in this system. 

Lastly let us discuss the equation of the level function associated with the present curvature flow 
equation. This method was used first by Ohta et a/[33|] in a special problem of the ordering process 



in a NCOP system with h = and has been developed by mathematicians independently: Let 
{u(r,t)} be the fictitious scalar field, the interface being defined by u(r,t) = 0. Then the interface 
normal velocity is represented as v n = |Vn|~ 1 ((?M/9t) u= o. Suppose we have an explicit curvature 
flow equation, 

v n (a)=F({H(a')}), (64) 
where T may be a functional. Here the mean curvature is given by 

H — —V • n = -OVwI^Vs • V< =0 , (65) 

where Vs = V — n(n -V). If we extend these equations to the bulk (u ^ 0), that makes a closure 
of Mr,*)}. 

Thus the level function equation for our case is written as 
du D£ . 2 _ 2 w IV7 ,_i 



dt 3A(t) 



Vu| y(-A(t)'Vs i )[|VM|- 1 V s • Vu + H}. (66) 
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For a small value of u we have a formula, 



BO 



H(u) - H(0) ~ u [\Vu\-\H 2 - 2K) + V s 2 | Vu\~ l } u=Q) (67) 

where .K" is the Gauss curvature defined by K — 1/R 1 R 2 using the radii of the principal curvature, 
Ri and R 2 . Then H in the right hand side of Eq. (|66"D may be replaced by 



H = H(0) + u - 2K)] u=0 . (68) 

The second term in the right hand side of Eq.([n]) canceled out after the surface integration. Note 
that the quantity H 2 — 2K (= R\~ 2 + -R2 2 ) is positive definite. Especially for a homogeneous 
bicontinuous phase of a symmetric mixture, this can be estimated as 

H 2 - 2K ~ -2K ~ 2A(t) 2 , (69) 

or alternatively we may use it as a definition of the parameter X(t). 

Thus the level function equation for a COP system has very tough form and has several 
difficulties|3^] to be solved. None of satisfactory analyses on it has been found so far.pBl Bll [M] 
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A on the pseudo-conformality 

Let f{u\,U2) be an arbitrary continuous function on S and suppose a harmonic function satisfying 
the Laplace equation 

(V s 2 + —a-e 3 — U( Ul ,u 2 ,u 3 ) = 0, (70) 
v e 3 du 3 du 3 ' 

with the boundary condition 

tp(ui,u 2 , 0) = f(ui,u 2 ). (71) 

Let us generalize the reference frame on S defined in §3 to the bulk, i.e., to the region 113 7^ in 
the same manner. Then by using the Green function method a formal solution for the harmonic 
function if> is given by 

Wtll ,„ 2 , M3 ) = - j f_|^ o/ ( ')rfi', (72) 

where 

G(u[, u' 2 , u' 3 ; «i, u 2 , us) = G + - G~, (73) 

and 

G = MR - «0 2 + ¥2 - U2) 2 + K t us) 2 } 1 ' 2 ' (74) 
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which satisfy a Poisson equation 

d 2 d 2 d 2 

W + du~i + d^ d± = ~ 6 ^' 1 ~ iLl)6 ^' 2 ~ T %) ' (75) 

Note that the reference frame is defined on taking the origin at the point just under the probe charge 
in §3. However, once the new Euclidean coordinate system is defined, one may shift the origin to 
any point. Thus, the reference coordinate (ui,u 2 ,u 3 ) is also used for the point of the probe charge 
in the above equations. 

Now combining Eqs.fl70D,fl72|) and (75 ) it is straightforward to find 



(£ + = (76) 

in the limit, u 3 — > 0. Here the following relation for small u 3 , 

1 d dG ± TT dG ± d 2 G ± d 2 G ± 

— ^— es^ — = H— — + — — ~ -^-g-, 77 

e 3 ou 3 ou 3 ou 3 du 3 du 3 

is used. Of course this approximation is exact for a minimal surface (H = 0). In general one may 
use the fact that a = —(dG ± /du' 3 ) u ' 3= o, which has an almost 5-function shape, satisfies 

da 1 _ 8 2 a 1 
ou 3 u 3 du 3 u 3 2 

and so on, when u 3 is sufficiently small. 

Using the same argument repeatedly on the functions Vs 2 /, Vs 4 / we obtain 

for an arbitrary integer n. This is the result desired in §3. 

It should be noted that this does not mean a conformal transformation on S, because the 
definition of the reference coordinate itself depends on the probed point. Further the reference 
coordinate u^) itself is not an orthogonal system on S in general. Then, strictly speaking, the 
operator d 2 /du 2 + d 2 jdu^ is not an invariant Laplacian but merely a Taylor expansion operator 
around the referred point. 

B a solvable example in d = 2 

It is instructive to discuss a solvable problem in the present scheme, though it is stupid practically. 
Let us consider a circle of radius R and assume a probe charge is located at a distance z from it. A 
new curvilinear coordinate {ux^Uz) may be defined by 

ui = (R + z)(p, u 2 = r-R, (80) 

where (r, ip) is the usual polar coordinate. The two dimensional Laplacian is written in this coor- 
dinate as 

~l/dlddd\ . . 

V = -(« ^ — + fi-Hiir-), 81 

gi K ou\giOUi 0U2 0U2' 



15 



where 

9l = = , gj, = iVtial" 1 = 1. (82) 

Let us derive the refraction law of the flux D in the Euclidean space (ui, u 2 ) due to the anisotropic 
dielectric tensor (ei, e 2 ) = (gi^ 1 , gi), which depends on u 2 only. Let a be the angle between the flux 
line and the wi-axis. Then we find that 

ei tan a = constant, (83) 
because of the continuity relations across a boundary line parallel to the Mi-axis, 

n . ^, . , -Dcosct D'cosa' , nA . 
Dsina = D sin a , = , (84) 

for the perpendicular and the parallel components, respectively. By using Eq. (|83|) the flux equation 
is given by 

du 2 tana R + u 2l 

— = = — — - — tana , (85) 

aui t\ H + z 

or 

log(J2 + u 2 ) = \og(R + z)- f^V (86) 
where «o is the initial angle, which defines the reference coordinate by 

ui = zj tana . (87) 
On putting u 2 = in Eq . (^6|) we find the mapping, U\ — > U\, i.e., 

Sl = (R + z)\o Z g (l + z/Rf l = \o g (l I z/Rf- (88) 

Of course this is a multi-valued mapping, (— ir, tt) — > (—00, 00). Then the induced charge density is 
given by 

1 z x 00 1 

where a new parameter x = \og(R + z) is used. By using some formulae this yields 

1 sinhx 

CindCfy? = 7 dip 

In cosh x — COS if 

1 {2R + z)z 

2^R 2 + (R + zf - 2R{R + z) cos <p V ' (90) 

That is the well-known result obtained directly with use of a mirror image method. 

Thus the multi- valued mapping has no problem in the present reference frame method. However, 
the cut-off procedure causes literally the cut-off of the tail part of the induced charge density. Here 
it is wrong to conclude that this effect would become infinitesimal in the limit z — > because the 
cut-off charge is of order of z/X(t). Note that the final formula for the equivalent simple layer is 
obtained after differentiating with respect to z. The weight for the contribution from the remote 
points becomes ~ du/u d , lacking the factor z. That estimates ~ 1/A(£) for the cut-off effect. Thus 
the cut-off loss must be recovered by some trick, such as the mean field term H/X(t) in Eq.( 
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